home *** CD-ROM | disk | FTP | other *** search
/ CD ROM Paradise Collection 4 / CD ROM Paradise Collection 4 1995 Nov.iso / hamradio / sgp4_pl2.zip / SGP4TST2.PAS < prev    next >
Pascal/Delphi Source File  |  1992-09-03  |  2KB  |  48 lines
  1. Program SGP4_Test;
  2. {$N+}
  3.   Uses CRT,
  4.        SGP_Init,SGP_Conv,
  5.        SGP_Math,SGP_Time,
  6.        SGP4SDP4;
  7.  
  8. var
  9.   satnumber,interval,i    : integer;
  10.   delta,time,tsince,k1,k2 : double;
  11.   pos,vel                 : vector;
  12.  
  13. BEGIN
  14.  
  15.   sat_data[1,1] := '1 88888U          80275.98708465  .00073094  13844-3  66816-4 0    8 ';
  16.   sat_data[1,2] := '2 88888  72.8435 115.9689 0086731  52.6988 110.5714 16.05824518  105 ';
  17.   sat_data[2,1] := '1 11801U          80230.29629788  .01431103  00000-0  14311-1        ';
  18.   sat_data[2,2] := '2 11801  46.7916 230.4354 7318036  47.4722  10.4117  2.28537848      ';
  19.  
  20.   delta := 360;
  21.   for satnumber := 1 to 2 do
  22.     begin
  23.     ClrScr;
  24.     Writeln(sat_data[satnumber,1]);
  25.     Writeln(sat_data[satnumber,2]);
  26.     Writeln;
  27.     Writeln('     TSINCE              X                Y                Z');
  28.     GotoXY(1,12);
  29.     Writeln('                      XDOT             YDOT             ZDOT');
  30.     Convert_Satellite_Data(satnumber);
  31.     time := Julian_Date_of_Epoch(epoch);
  32.     tsince := 0;
  33.     for interval := 0 to 4 do
  34.       begin
  35.       SGP(time,pos,vel);
  36.       Convert_Sat_State(pos,vel);
  37.       GotoXY(1,6+interval);
  38.       Writeln(tsince:16:8,pos[1]:17:8,pos[2]:17:8,pos[3]:17:8);
  39.       GotoXY(1,14+interval);
  40.       Writeln('                ',vel[1]:17:8,vel[2]:17:8,vel[3]:17:8);
  41.       tsince := tsince + delta;
  42.       time := time + delta/xmnpda;
  43.       end; {for int}
  44.     repeat until keypressed;
  45.     end; {for satnumber}
  46.  
  47. END.
  48.